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Abstract 

Magnetization, resistivity and specific heat measurements were performed on solution-grown, 



i. 

c/2 . single crystals of RFe2Zn2o and RCo2Zn2o(R = Y, Nd, Sm, Gd - Lu). Whereas LuCo2Zn2o and 



YCo2Zn2o manifest unremarkable, metallic behavior, LuFe2Zn2o and YFe2Zn2o manifest behaviors 
characteristic of nearly ferromagnetic Fermi liquids. When the well defined 4/ local moments 
(Gd^^ - Tm^+) were embedded into this strongly polarizable host, they all manifest enhanced, 
ferromagnetic ordering and the values of Tc for RFe2Zn2o (R = Gd - Tm) scale with the de Gennes 
factor. In addition, the data on the RFe2Zn2o compounds indicate a small crystal electronic field 
(CEF) effect, compare with the interaction energy scale. On the other hand, the local moment 
bearing members of RCo2Zn2o (R = Nd, Sm, Gd - Tm) manifest weak magnetic interactions and 



I the magnetic properties for R = Dy - Tm members are strongly influenced by the CEF effect on the 



R ions. The magnetic anisotropy and specific heat data for the Co series were used to determine 
the CEF coefficient of R ion in the cubic coordination. These CEF coefficients, determined for 
the Co series, are consistent with the magnetic anisotropy and specific heat data for the Fe series, 
which indicates similar CEF effects for the Fe and Co series. Such analysis, combined with specific 
heat and resistivity data, indicates that for R = Tb - Ho, the CEF splitting scale is smaller than 
their Tc values, whereas for ErFe2Zn2o and TmFe2Zn2o the 4/ electron loses part of its full Hund's 
rule ground state degeneracy above Tq- YbFe2Zn2o and YbCo2Zn2o manifest typical, but distinct 
heavy Fermion behaviors associated with different Kondo temperatures. 

PACS numbers: 75.30.Gw, 75.50.Ee, 75.50.Cc 
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I. INTRODUCTION 



Intermetallic compounds consisting of rare earth and transition metals, as well as metal- 
loids, have versatile magnetic properties. Compounds with itinerant d electrons are of 
particular interest when they are in the vicinity the Stoner transition: such systems, charac- 
terized as nearly or weakly ferromagnet, manifest strongly correlated electronic properties.!^ 
On the other hand, heavy rare earth ions manifest magnetic versatility associated with the 
4/ electrons: null magnetism (Y^"*" or Lu'^^), pure spin, local moment magnetism (Gd^"*"), 
potentially anisotropic, crystal electric field (CEF) split, local moment magnetism (Tb^+ - 
Tm^+), and more exotic magnetism: Yb ions may hybridize with conduction electrons and 
manifest so-called heavy fermion behavior. Needless to say, series of materials that combine 
these interesting versatilities have attracted the attention of physicists. For example, the 
binary RC02 (R = rare earth) compounds, with the nearly ferromagnetic (FM) end members 
YC02 and LuCo2, and the local moment, FM members (R = Pr, Nd, Gd - Tm), have been 
studied for more than 35 years. 



Discovered by Nasch et. al 



a 

6!], the RT2Zn2o (R = rare earth, T = transition metal in 
ihe Fe, Co and Ni groups) series of compounds crystallize in the cubic CeCr2Al2o structure 
3, y, yl- The R and T ions occupy their own single, unique, crystallographic sites with 
cubic and trigonal point symmetry respectively, whereas the Zn ions have three unique 
crystallographic sites (Fig. [T]a). Both the R and the T ions are fully surrounded by the 
shells of Zn that consist of the nearest neighbors (NNs) and the next nearest neighbors 
(NNNs). This means that there are no R-R, T-T or R-T nearest neighbor and the shortest 
R-R spacing is ~ 6 A. For the R ions, a Frank-Kasper polyhedron with coordinate number 
(CN) 16 is formed from 4 Zn NNs on the 16c site and 12 Zn NNNs on the 96^' site (Fig. [1] 
b). 

Recently, studies of the thermodynamic and transport properties of these intermetallics 
have revealed varied magnetic behavior. jiol. 11 . 12 . 13 | YFe2Zn2o and LuFe2Zn2o are 



archetypical examples of nearly ferromagnetic Fermi liquids (NFFL) with Stoner enhance- 
ment factors oi Z = 0.88 (where Xt=o = XPauu/i.^ — Z)). By embedding large, Heisenberg 
type moments associated with Gd^^ ions in this highly polarizable medium, GdFe2Zn2o man- 
ifests highly enhanced FM order. On the other hand, GdCo2Zn2o manifests ordinary, low 
temperature, antiferromagnetic (AFM) order, correspondent to the 'normal metal' behavior 
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FIG. 1: (Color online) (a) The cubic unit cell of RT2Zn2o- (b) The CN16 Frank-Kasper polyhedron 
of rare earth ions. 

of the conduction electron host, YCo2Zn2o-|l^ In addition to these interesting, 3d electron 
and local moment properties, six related YbT2Zn2o compounds (T = Fe, Co, Ru, Rh, Os 
and Ir) show heavy fermion ground states, associated with different Kondo temperatures 



(Tk) and Yb ion degeneracies. [11] 



Given the similarities and differences between the RFe2Zn2o and RCo2Zn2o (R = Gd, Y, 
Lu) series, it becomes important to study all of the R = Y, Gd - Lu members in detail. 
A comparative study of the RFe2Zn2o and RCo2Zn2o series will help to further understand 
the magnetic interaction between the local moments by means of the strongly polarizable 
medium, particularly with the crystal electronic field (CEF) effect associated with non- 
zero orbital angular momentum. Furthermore, given the very similar CN-16 Frank-Kasper 
polyhedron for R ions, as well as the less than 2% difference of lattice constants for the 
whole RT2Zn2o families, the study of the CEF effect on these local moment members will 
also help in the understanding of the varied heavy fermion states of YbT2Zn2o, which were 
thought to be due to the competition between temperature scales associated with the CEF 



splitting and the Kondo effect. |lll| 



In this paper, we present the results of magnetization, heat capacity and resistivity mea- 
surements on RFe2Zn2o and RCo2Zn2o (R = Y, Nd, Sm, Gd - Lu) compounds. Com- 
pared with the weakly correlated behaviors for YCo2Zn2o and LuCo2Zn2o, YFe2Zn2o and 
LuFe2Zn2o, manifest clear, NFFL behaviors associated with the spin fluctuation of the itin- 
erant electrons. For the RFe2Zn2o compounds (R = Gd - Tm), the well-defined, local moment 



members all manifest enhanced FM ordering with Tq values that roughly scale with the de 
Gennes factor. Their anomalous, temperature dependent susceptibility and resistivity can 
be explained as the result of local moments embedded in a NFFL host. In contrast, for the 
RCo2Zn2o series, only Gd and Tb members manifest AFM ordering above 2 K, and the mag- 
netic properties for R = Dy - Tm clearly manifest features associated with single ion CEF 
effects on the R ions in the cubic symmetry coordination. For the R = Tb - Tm members 
in the Co series, the CEF parameters can be determined from the magnetic anisotropy and 
the specific heat data, and are roughly consistent with the results of calculations using a 
point charge model. For the Fe series, the R = Tb - Tm members show moderate magnetic 
anisotropy in their ordered states, mainly due to the CEF effect on the R ions, which is 
consistent with the magnetic anisotropy found for the Co members. These results, as well 
as the analysis of the heat capacity and resistivity data, indicate that the FM state emerges 
from the fully degenerate Hund's rule ground state for RFe2Zn2o (R = Gd - Ho), whereas 
ErFe2Zn2o and TmFe2Zn2o manifests CEF splitting above their Curie temperatures. 



II. EXPERIMENTAL METHODS 



Single crystals of RFe2Zn2o and RCo2Zn2o (R = Y, Gd - Lu) were grown from high 
temperature, Zn rich solutions with the initial concentration of starting elements being 

nn 

R:(Fe/Co):Zn = 2:4:94, as described previously |10L Il4|. The small amount of residual Zn 
on the resulting crystals was removed by submerging the crystals into an ultrasonic bath of 
0.5 vol. % HCl in H2O for 0.5 - 1 hour. Attempts were made to grow even lighter rare earth 
members (R = Eu, Nd and Ce) for the Fe series, but no ternary compound was found to 
form. On the other hand, for the Co series, single crystals for the R = Nd and Sm members 
were grown using the same growth method. Room temperature, powder X-ray diffraction 
measurements, using Cu Ka radiation, were performed on the samples with additional Si 
powder (a = 5.43088 A) used as standard. Figure [2] shows the lattice constants for the two 
series of compounds, obtained by using the Rietica, Rietveld refinement program, plotted 
with respect to the effective radius of R^"*" with CN = 9 [15], since the data is absent for larger 
CN. The variation of the lattice constant manifests the well-known, lanthanide contraction 
for R = Gd - Lu with no evident deviation for R = Yb. However, the relatively larger 
lattice constants for YFe2Zn2o and YCo2Zn2o indicate that, with this large CN, the effective 
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FIG. 2: (Color online) The lattice constants (a) for RFe2Zn2o and RCo2Zn2o versus the radius of 
the trivalent rare earth ion with CN = 9. [iq] The error bars were estimated from the standard 
variation of four seperate measurement made on one batch of sample. 



ionic radii of is between Tb'^"'' and Dy'^+ rather than between Dy^+ and Ho'^+. This 
deviation for ions is not unprecedented in the isostructural compounds RRu2Zn2o 
and RMn2ln3;Zn2o-x ^^^^ similar structure compounds RC02 17]. Additional 

single crystal X-ray diffraction measurements were preformed on R = Gd, Tb, Er and Lu 
members of the RFe2Zn2o series and demonstrated full occupancy on all cr y st allogr aphic 



sites (within the detection errors) and the same lattice as the powder X-ray values. [l8j This 
result is complicated by the difficulty in resolving the difference between Fe and Zn elements 
with very similar X-ray scattering strengths. 

It is worth noting though, that for the local moment bearing members in Fe series, single 
crystals obtained from different ratios of starting element concentrations manifest detectably 
different magnetic ordering temperatures. These differences, tentatively associated with 
very subtle variations of element occupancy on the crystallographic sites[13], are related to 
the extreme sensitivity to the small disorder for compounds with such strongly correlated 
conduction electron backgrounds. A detail discussion of this is presented in Appendix A. 

Measurements of electrical resistivity were preformed by using a standard AC, four-probe 
technique. The samples were cut into bars, with typical lengths of 2-3 mm, along their 
crystallographic [110] direction. These bars were measured in a Quantum Design physical 
properties measurement system, PPMS-14 and/or PPMS-9 instruments (T = 1.85 - 310 
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K) with / = 16 Hz, 7 = 3- 0.3 mA. Temperature dependent specific heat measurements 
were also performed in these PPMS units, using the heat capacity option. In order to 
study the magnetic part of specific heat, the contribution from the conduction electrons and 
lattice vibrations was subtracted. This was performed by subtracting off the specific heat of 
the non-magnetic members YT2Zn2o or LuT2Zn2o [Cmag = Cp(RT2Zn2o) — Cp(Lu/YT2Zn2o)]. 
The magnetic part of entropy can then be estimated as Sm = J '~^p^^^^'^^^o^~^^^^/^^^'^^^o^ (YJ"^ 
in which the specific heat data below 1.8 K (0.4 K for TmCo2Zn2o) were approximated by 
a linear extrapolation to Cp = at K. 

Magnetization measurements were carried out in Quantum Design magnetic properties 
measurement system (MPMS), superconducting quantum interface device (SQUID) magne- 
tometers in varied applied fields { H < 55 kOe) and temperatures (T = 1.85 - 375 K). In 
the measurements of magnetization for the FM-ordered, RFe2Zn2o compounds, the effects 
of the dernagnetizing field are relatively small, because of the dilute nature of the magnetic 
moments: |l2 . the maximum demagnetizing field is about 3000 Oe for a plate-like shaped 
sample. Even so, in order to obtain accurate magnetization isotherms near Tq, the samples 
were cut and measured with the applied field along their long axis so as to minimize this 
already small demagnetization effect. 



III. RESULTS 



We start characterizing the compounds with the non-magnetic rare earth ions of the 
series: Y(Lu)Fe2Zn2o and Y(Lu)Co2Zn2o. Without any 4/ electronic magnetism, these 
compounds manifest the electronic and magnetic properties associated with the conduction 
electron background of each series. Next, we will introduce the two series of compounds 
with well-defined 4/ local moments: R = Gd - Tm. We will examine the magnetization and 
specific heat data for the Co series at first. Then an overview of the magnetic properties 
for the Fe series will be presented next. After that, the magnetization, specific heat and 
resistivity data will be presented for each Fe member separately. Finally, similar data for 
the R = Yb heavy fermion compounds, YbFe2Zn2o and YbCo2Zn2o will be presented. 



6 




I 1 I 1 I n ^^wel I , I 

100 200 300 400 



T(K) 

FIG. 3: (Color online) (a) Temperature dependent magnetization M divided by the applied field 
H for YFe2Zn2o and LuFe2Zn2o as well as their Co analogues for i7 = 10 kOe and 50 kOe. Inset: a 
blow-up plot at low temperature, (b) H/M for YFe2Zn2o and LuFe2Zn2o. The solid lines present 
the modified Curie- Weiss [xiT] = C/{T — 6c) + Xo] fit for the data above 100 K. Inset: field 
dependent magnetization at 1.85 K. 



A. Y(Lu)Fe2Zn2o and Y(Lu)Co2Zn2o 

Temperature dependent magnetization data (divided by the applied field) for 
Y(Lu)Fe2Zn2o and Y(Lu)Co2Zn2o are shown in Fig. [3] a. The Fe members manifest similar, 
strongly enhanced, temperature-dependent paramagnetic signals, whereas the Co members 
manifest essentially temperature- independent, Pauli paramagnetic signals. The low temper- 
ature features for the two Fe compounds are shown in the inset of Fig. [31 In the applied 
field of 10 kOe, the magnetization signals of both Fe members show a faint maximum below 
10 K, whereas the high magnetic field (50 kOe) suppresses the lowest temperature M/H 
values, as well as the maximum. In our experience on the measurements of different batches 
of samples, these low temperature features are moderately sample-dependent (different sam- 
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pies may show 20% different magnetization signal and 2-3 K difference in the temperature 
of the maximum, Tmax)- Nevertheless, the maximum of the temperature dependent suscep- 
tibility, x(T), is a common feature in the NFFLs, such as Pd [20], YC02 and LuCo2 [21], as 
221 ]. although quantitative calculation of x{T) still presents a challenge even 



well as TiBe2 



for the simplest case of Pd 23|, |2J]. The field suppression of the magnetization (and T„ 



at low temperature is not attributed to the possible existence of a paramagnetic impurity 
contribution (which would contribute more to the value of M/H at lower temperature and 
lower field, and therefore suppress the maximum of M/H in lower field), but, as discussed 
below, to the intrinsic variation of x = dM/ dH with respect to H at different temperatures. 

Figure Ob shows that above a characteristic temperature (T* ~ 50 K), the susceptibil- 
ity of YFe2Zn2o and LuFe2Zn2o can be approximately fitted by a Curie- Weiss (CW) term 
[x(T) = C/(T — 9c)] plus a temperature- independent term (xo)- The values of effective 
moment (fieff), 9c and are extracted as 1.0 yU^/Fe, -16 3.8 x 10~^emu/mol and 1.1 
/ifi/Fe, -33 3.4 x 10~^emu/mol for YFe2Zn2o and LuFe2Zn2o, respectively. These values of 
IXeff are significantly larger than the estimated induced moment of Fe site in the FM ground 



state of GdFe2Zn2o, ~ 0.35uR/Fe. [lOj Such apparent CW-like behavior was a! 



in other NFFL systems. 
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251 ] In the context of the spin fluctuation model \ 



so observed 
3|, itinerant 



electronic systems can manifest CW-like behavior with a Curie constant related to the local 
amplitude of the spin fluctuation. The magnetization data at the base temperature (1.85 
K) show nearly linear dependent with the applied field (Inset in Fig. [3]b), which is distinct 
from the Brillouin function type of magnetization curves associated with local moments. 

In order to better understand the variation of the maximum in temperature dependent 
M/H data for YFe2Zn2o and LuFe2Zn2o, M{T) and M{H) measurements were performed 
on LuFe2Zn2o for varied applied field and temperature respectively. Figure H] shows that 
the magnetic field suppresses the values of M/H, as well as the maximum of M/H, which 
disappears when H > 20 kOe. The values of the field dependent susceptibility, x{H), were 
extracted as x{H) = ^^^^l^j'^^^^^ with AH = 1000 Oe from the M{H) data at varied 
temperature (Figure [5]). For T > 10 K, the values of x{H) monotonically decrease with 
increase H, whereas a local maximum appears around 20 kOe in the data sets as T < 7 K. 
This critical temperature (~ 7 K) is close to the T^ax] the maximum of xi^) {H = 20 kOe) 
is also close to the suppression field determined by Fig. HI This curious, field dependent, 
susceptibility at varied temperature is reminiscent to the one of TiBe2, albeit the amplitude 
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FIG. 4: M/H for LuFe2Zn2o. From right to left: H = 2 kOe, 5 kOe, 10 kOe, 20 kOe and 30 kOe. 
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FIG. 5: x{H) for LuFe2Zn2o at varied temperature. 
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261] In the case of TiBe2, the reason of 



of local maximum in x{H) is much smaller, 
anomalous field- dependent magnetization is also still not clear. [2j 

Figure M presents the low temperature specific heat data for YFe2Zn2o and LuFe2Zn2o, 
as well, as for the Co analogues, plotted as Cp/T versus T^. All four compounds manifest 
clear Fermi liquid behavior {Cp = 7T + (3T^). The similar (3 values (represented as the 
slopes of the data sets in the plot, ~ 1.2mJ/molK^) indicate the similar Debye temperatures 
for these 4 compounds (~ 340 K [13]), consistent with their similar molar mass, similar 
composition and similar lattice parameters. On the other hand, the over 2.5 times larger 
values of electronic specific heat (7) of the Fe members indicate a larger density of states 
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FIG. 6: (Color online) Low temperature specific heat data of YFe2Zn2o and LuFe2Zn2o (plotted as 
Cp/T versus T^), as well as the Co analogues. 



at Fermi level [Nj Ef)] , compared to the Co analogues (consistent with the band structure 
calculation results |l3[|). 

The values of the electronic specific heat and magnetic susceptibility can be employed to 
estimate the Stoner enhancement factor, Z, in the context of the Stoner theory: that is, the 
static susceptibility is enhanced by , whereas the electronic specific heat is not. [lo| 
The estimated Z values of YFe2Zn2o and LuFe2Zn2o are 0.88, 0.89, respectively, comparable 
with the estimated values of the canonical NFFL systems: Pd: 0.83, and YC02: 0.75. |28| 

The temperature dependent electrical resistivity data for YFe2Zn2o and LuFe2Zn2o are 
larger than that for the Co analogues over the whole temperature range (Fig. [7] a). This is 
not unexpected for a NFFL since the spin fluctuations will affect the scattering process of the 
conduction electrons, which leads to an additional contribution to the resistivity. In order 
to study the spin fluctuation contribution to the resistivity, the total electrical resistivity 
p(T) is assumed to be: 

p{T) = p^ + Ppn{T) + psf{T), (1) 

where the first, second and third terms represent residual, phonon and spin fluctuation 
reisitivity, respectively. Assuming the phonon scattering contribution, ppfi{T), is essentially 
same for the Fe and Co analogues (as suggested by the similar (3 terms), then, the spin 
fluctuation scattering contribution, psf{T), can be estimated as: 
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FIG. 7: (Color online) (a): Temperature dependent resistivity of YFe2Zn2o and LuFe2Zn2o, as well 
as their Co analogues, (b): estimated spin fluctuation contribution to the resistivity for YFe2Zn2o 
and LuFe2Zn2o. The error bars were estimated as ±10% of the values of the resistivity for YCo2Zn2o 
and LuCo2Zn2o respectively. 



Psf{T) — (p - Po)Y/LuFe2Zn2o " (P " Po) Y/LuCoaZnao • (2) 

Shown in Fig. [7)d, Psf{T) for these two compounds increases with temperature and is 
close to a saturated value (10 fiQcm) at 300 K, within the accuracy of the measurements. 

The analysis of the low temperature resistivity data reveals a quadratic, standard Fermi 
liquid, behavior [p(T) = po + AT"^] for all 4 compounds (Fig. [8]). The A values of the Fe 
compounds are about 7 times larger than the two Co analogues. This result is consistent 
with the 2.5 times larger 7 values of the Fe compounds, in the context of the Fermi liquid 
theory, meaning A is proportional to the square of the effective mass of the quasi-particles 
due to the strong correlation effect, whereas 7 is proportional to the effective mass. In the 
point of view of spin fluctuation theory, nearly FM metals manifest Fermi liquid behaviors 
at low temperature region with A values enhanced by spin fluctuations.!^ 
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FIG. 8: (Color online) p versus for YFe2Zn2o and LuFe2Zn2o, as well as their Co analogues. 
The solid lines present the linear fit of the data sets from 2 K to 9 K. 
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FIG. 9: (Color online) Temperature dependent magnetization of RCo2Zn2o (R = Nd and Sm) 
compounds, divided by applied field H = 10000 Oe. Inset: applied field {H = 10000 Oe) divided 
by the magnetizations of RCo2Zn2o (R = Nd and Sm) as a function of temperature. 

B. RCo2Zn2o (R = Nd, Sm, Gd - Tm) 

Before discussing the heavy rare earth compounds (R = Gd - Yb), the results of ther- 
modynamic measurement on NdCo2Zn2o and SmCo2Zn2o will be briefly examined. Figure [9] 
shows the temperature dependent magnetization data (divided by the applied field H = 1000 
Oe) for NdCo2Zn2o and SmCo2Zn2o- Neither of them manifest any sign of magnetic order- 
ing above 2 K. The temperature dependent H/M for NdCo2Zn2o shows a CW behavior 
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FIG. 10: (Color online) Temperature dependent specific heat for RCo2Zn2o (R = Nd, Sm and Y). 



[x{T) = C/{T - Oc) + Xo] with ^^ff = 3.7/iB, Be = -2.3 K and Xo = 6.8 x 10"^ emu/moh 
The value of the effective moment is close to the theoretical values for the Hund's rule ground 
state of the 4/ electrons of Nd^"*" ion (3.6/iB). On the other hand, though the magnetization 
of SmCo2Zn2o drops with increase temperature, it does so in a distinctly non-CW manner. 
This behavior is not unexpected in Sm containing compounds 29(], and is most likely due 
to the thermal population of the first excited Hund's rule multiplet. 

Specific heat data for NdCo2Zn2o and SmCo2Zn2o are shown in Fig. [10] along with data 
for YCo2Zn2o for comparison. The low temperature upturn in the NdCo2Zn2o data below 2 
K may be due to a lower temperature magnetic ordering or a Schottky anomaly due to the 
CEF splitting. The specific heat data for SmCo2Zn2o manifest a broad peak around 4 K, 
which is most likely due to the CEF splitting of the Hund's rule ground state of Sm^"*". Both 
NdCo2Zn2o and SmCo2Zn2o data increase much faster above 10 K, and remain more than 10 
J/mol K larger above 25 K (not shown here), compared with the data for the non-magnetic 
analogue YCo2Zn2o. On the other hand, the calculated results of the CEF splitting for the 
Hund's rule ground state of Nd^^ ion in a point charge model show the splitting energy 
levels within 25 K (see Table HTl below). This large difference indicates that, at this point, 
the magnetic part of Cp for NdCo2Zn2o and SmCo2Zn2o cannot be well estimated, since the 
Cp data of YCo2Zn2o is not an adequate approximation of the non-magnetic background 
and unfortunately the LaCo2Zn2o analogue did not form. 

Temperature dependent magnetization data (divided by the applied field H = 1000 Oe) 
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FIG. 11: (Color online) Temperature dependent magnetization of RCo2Zn2o (R = Gd - Yb) com- 
pounds, divided by applied field H = 1000 Oe. 



for RCo2Zn2o (R = Gd - Yb) are presented in Fig. [TTl In addition to the previously reported, 
AFM ordered GdCo2Zn2o with the Neel temperature Tn = 5.7 ± 0.1 kQ, Q, TbCo2Zn2o 
also shows AFM ordering with Tn = 2.5 ± 0.1 K, which also clearly manifests itself in the 
specific heat data (shown below in Fig. [T3l) . The rest of the members (R = Dy - Yb) do not 
show magnetic ordering above 2 K. Due to the relatively low density of state at Fermi level 



[N{Ef)] for the Y and Lu analogues |10l. Il3| and large R-R separation, such low temperature 
magnetic ordering for the 4/ local moments coupled via the Ruderman-Kitter-Kasaya-Yosida 
(RKKY) interaction is not unexpected. 

Figure fT2] shows the temperature dependent H/M for R = Gd - Tm and Yb members 
of the RCo2Zn2o series. All the members, including YbCo2Zn2o, manifest CW behavior 
[x(T) = C/(T — 6c) + Xo] with negligible small Xo 2 x lO^^emu/mol) and the values 
of /ie// close to the theoretical values for the Hund's ground state of the 4/ electronic 
configurations; all the values of 6c are close to 0, consistent with the low magnetic ordering 
temperatures (Table [T]). 

The specific heat data for RCo2Zn2o (R = Gd - Tm , Y and Lu), as well as the pseudo- 
ternary compound Tbo.5Yo.5Go2Zn2o are presented in Fig. [131 In addition to the previously 
studied GdCo2Zn2o, the specific heat data for TbCo2Zn2o manifests a A-type of anomaly 
with a peak position at 2.5 K, the AFM ordering temperature. In addition to this peak, 
the Cp data also show a broad shoulder above 2.5 K, which is due to the CEF splitting 
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FIG. 12: (Color online) Applied field {H = 10000 Oe) divided by the magnetizations of RCo2Zn2o 
(R = Gd - Yb) as a function of temperature. 



TABLE I: Paramagnetic Curie temperature, 9c (with ±0.1 K errors) and effective moment, ^e// 
[from the CW fit of x{T) from 50 K to 300 K]; Neel temperature, Tn for RCo2Zn2o compounds (R 
= Nd, Gd - Yb). 
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above the magnetic ordering temperature. This anomaly, associated with CEF splitting 
of the 4/ electronic configuration of Tb^"*", manifests itself more clearly in the Cp data for 
Tbo.5Yo.5Co2Zn2o: when Tn is suppressed to well below 2 K, the Cp data show a Schottky 
anomaly with a peak position ~ 3 K. The magnetic part of entropy for TbCo2Zn2o is shown 
in the inset to Fig. [131 Approximately 50 % of the total magnetic entropy is recovered by 
Tn, and by 15 K the full Sm = -Rlnl3 is recovered {R is gas constant). This is consistent 
with a very small, total CEF splitting in the these compounds, associated with the highly 
symmetric environment of the R ions. For the rest of the members, R = Dy - Tm, the 
specific data show broad, Schottky-type of anomaly below 10 K, as shown in the insets of 
Fig. [16], [T71 [18] and [19] (shown below). The low temperature upturn for DyCo2Zn2o data 
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FIG. 13: (Color online) Temperature dependent specific heat for RCo2Zn2o (R = Gd - Tm, Y and 
Lu), as well as Tbo.5Yo.5Co2Zn2o- Inset: temperature dependent magnetic entropy for TbCo2Zn2o- 
The dashed line presents the entropy of the full Hund's ground state of Tb"*"^. 
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FIG. 14: (Color online) Normalized magnetic part of entropy for RCo2Zn2o (R = Dy - Tm) as well 
as for Tbo.5Yo.5Co2Zn2o (in units of per mole R^+). The error bars were estimated from the ±1 % 
of the total entropy. 

below 2 K may indicate a magnetic ordering at lower temperature, whereas the upturn 
for TmCo2Zn2o data below 0.7 K may be associated with a magnetic ordering at very low 
temperature and/or a nuclear Schottky anomaly. 

The released, magnetic part of entropy above 2 K (above 0.4 K for the TmCo2Zn2o) 
are shown in Fig. [HI For R = Dy - Tm, there is an obvious deficit of magnetic entropy 
compared with the value associated with fully degenerated Hund's ground state, which 





16 



indicates unaccounted entropy below 2 K (0.4 K for TmCo2Zn2o) associated with low lying 
CEF levels and magnetic ordering. 

In order to better understand the magnetic properties for R = Tb - Tm members, the CEF 
effect acting on the R ions is evaluated by thermodynamic measurements. The single-ion 
Hamiltonian for the R'^^ is assumed to be the sum of the CEF term, an exchange interaction 
term and an external field term: 

Ti = HcEF + 'Hexc + 'Hext- (3) 

where Tiext = gjl^sJ- H, gj is Lande factor, J is the total angular momentum, and H is 
the external magnetic field. 

Since the rare earth ions are located in a cubic point symmetry, the CEF term, Hcef, 
can be written as: 

Hcef = B^Ol + 50^) + i?6°(0° - 2106^). (4) 

where operators are the well-known Stevens operators [soj, and and Bq are CEF 
parameters [31] . If one follows the work of Lea et al. [31], this expression can be written as: 

-Hcef = W[jr^{Ol + Wt) + ^-^{O^ - 2lOt)]. (5) 



where F4 and FQ are factors introduced by Lea et al. 311] and dependent with J, W is the 
energy scale, and x represents the relative importance of the 4th and 6th order terms. 

Noticing that the possible magnetic ordering temperatures are below 2 K for RCo2Zn2o 
(R = Dy - Tm), as well as for Tbo.5Yo.5Co2Zn2o, the exchange interaction term will be 
approximated as zero, an approximation that will be better for R = Tm than for R = Dy. 
Thus, the CEF parameters for different R ions were determined by fitting the magnetization 
at 2 K and the temperature dependent specific heat data. 

Figure [T514T9] show data and CEF fitting results of the magnetization at 2 K and the 
magnetic part of specific heat with the single ion Hamiltonian (ignoring the interaction 
term) for Tbo.5Yo.5Co2Zn2o and RCo2Zn2o (R = Dy - Tm). The specific heat data for all 
members are less than the one of YCo2Zn2o above 30 K, which is likely due to the errors 
associated with resolving the difference between the sample's total Cp and the relatively 
large nonmagnetic contribution. Therefore, the fittings of Cm were performed below 20 K. 
For R = Dy - Tm, the experimental magnetization data were slightly less than the calculated 
results. Such phenomena, more significant for R = Dy and Ho, are most likely due to the 
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FIG. 15: (Color online) Field dependent magnetization for Tbo.5Yo.5Co2Zn2o along three principle 
axes. The solid lines present the fitting results. Inset: magnetic part of specific heat. The solid 
and dashed line present the experimental and calculated result respectively. 
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FIG. 16: (Color online) Field dependent magnetization for DyCo2Zn2o along three principle axes. 
The solid lines present the fitting results. Inset: magnetic part of specific heat. The solid and 
dashed line present the experimental and calculated result respectively. 

still relevant AFM-type of interaction between the local moments. As shown in table [Tll the 
inferred W and x values for all 5 compounds are clustered in a narrow range: \W\ < 0.1, 
\x\ < 0.25. This result, indicating small energy scales of the CEF effect and relatively large 
Bq terms, are roughly consistent with the calculated results based on the point charge model 
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FIG. 17: (Color online) Field dependent magnetization for HoCo2Zn2o along three principle axes. 
The solid lines present the fitting results. Inset: magnetic part of specific heat. The solid and 
dashed line present the experimental and calculated result respectively. 




FIG. 18: (Color online) Field dependent magnetization for ErCo2Zn2o along three principle axes. 
The solid lines present the fitting results. Inset: magnetic part of specific heat. The solid and 
dashed line present the experimental and calculated result respectively. 



(see Appendix B). Furthermore, it should be noted that the signs of the Bq terms for the 
calculated results are all consistent with the experimental ones; this is not the case for the 
B2 terms. This behavior is not difficult to understand, as shown in the Appendix B, the 
contributions to the CEF splitting are mainly from the CN-16 Frank-Kasper polyhedron 
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FIG. 19: (Color online) Field dependent magnetization for TmCo2Zn2o along three principle axes. 
The solid lines present the fitting results. Inset: magnetic part of specific heat. The solid and 
dashed line present the experimental and calculated result respectively. 



TABLE II: Comparison of the CEF parameters of RCo2Zn2o compounds (R =Nd, Tb - Yb), 
determined from magnetization measurements to those calculated in a point charge model. 
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formed by 4 NN and 12 NNN Zn neighbors. For the term, the contributions cancel each 
other by the two sets of neighbors, whereas the contributions for the terms is the sum. 
Therefore, the Bq terms are relatively large and the calculated results are more reliable. 
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FIG. 20: (Color online) Temperature dependent magnetization of RFe2Zn2o (R = Gd - Tm), 
divided by applied field H = 1000 Oe. 



C. RFesZnao (R = Gd - Tm) 

Before discussing each of the compounds in this series separately, an overview of their 
temperature dependent magnetization data serves as a useful point of orientation. Figure 
[20] shows M/H versus T (the applied field H = 1000 Oe) for R = Gd - Tm members. In 
contrast to the Co series compounds, the Fe series compounds all manifest FM ground states 
with enhanced Tq values, which systematically decrease as R varies from Gd to Tm. Such 
enhanced FM ordering has been explained as the result of local moments embedded in the 
NFFL host, most clearly seen in YFe2Zn2o and LuFe2Zn2o- [10] This systematic variation 
of Tc on R is not unexpected for such heavy rare earth compounds when the magnetic 
interaction between the R ions are associated with the spin part of the Hund's ground state 
of 4/ electrons. 

The temperature dependent H/M data, approximately equaling inverse susceptibihties 
[l/x(r)] in the paramagnetic state, for R = Gd - Tm, as well as for YbFe2Zn2o, are shown 
in Fig. [211 Similar to GdFe2Zn2o \1Q\, the l/x(T) data sets for R = Tb - Tm follow the 
CW law [x(T) = C/(T — 6c) + Xo] at high temperatures, and deviate from the law when 
approaching their magnetic ordering. The effective moment (yUe//) and the paramagnetic 
Curie temperature {9c) data for these 6 compounds are listed in Table [TTTl All /ig// values 
are close to the theoretical value for the Hund's ground state of the trivalent 4/ electronic 
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FIG. 21: (Color online) Applied field (H = 1000 Oe) divided by the magnetizations of RFe2Zn2o 
(R = Gd - Tm) function of temperature. 

TABLE III: Residual resistivity ratio, RRR = R{300K) / R(2K); paramagnetic Curie temperature, 
6c (with ±0.5 K errors) and effective moment, /^e// (from the CW fit of xi^) from 100 K to 300 
K, except for GdFe2Zn2o, which was fitted from 200 K to 375 K Q); Curie temperature, Tc; and 
saturated moment at 55 kOe along the easy direction, ij,sat for RFe2Zn2o compounds (R = Gd - 
Yb). 
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configuration. 



1. TbFe2Zn2o 

Temperature dependent M/H, specific heat and resistivity data sets for TbFe2Zn2o are 
shown in Fig. [22l The M{T)/H data are consistent with FM order below 60 K, and the 
magnetic phase transition manifests itself as a faint feature in Cp data, indicating Tc = 
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FIG. 22: (a) Temperature dependent M/H for TbFe2Zn2o {H = 1000 Oe); (b) C^; (c) p and 
dp/dT. Upper inset : magnetic part of specific heat. Lower inset: magnetic entropy Sm- 

o 
E 

I 20 

-Si- 

5 10 





^300 



o 

E200 



o°ioo 



30 

"e 

" 20 

a 

<i 10 




1 1 ■ 1 1 1 1 

' yTbFe^Z^ a); 


O 

■ E 

3 5 
-(J 

■ 4 


1 

20 40 60 8 
T(K) > 


1 

25 
20 
E 15 
3 10 
5 

°C 


1 

20 40 60 80 
T(K1 


/ , 1 . 1 , 


1 



Q. 

■o 



20 



40 

T(K) 



60 



80 



0.4 
0.3 
0.2 ^ 
0.1 
.0.0 



56 ± 3 K. This is clearer in the magnetic part of specific heat {Cm) data (Fig. [22] upper 
inset) where Tq is taken as the position of the greateat upward slope. As shown in the lower 
inset to Fig. [22] b, at Tc, the magnetic entropy is close to the value for the full degeneracy 
of the Hund's ground state of Tb^"*", i?ln 13. As we shall see for the rest of the local moment 
members (R = Dy - Tm), the released magnetic entropy at Tc for the respective rare earth 
ion, is close to the full degeneracy value for their Hund's ground state, except for R = Tm. 
The p(T) data manifests a change in the slope, which could be seen even more clearly on in 
the dp/dT data, consistent with a Tc = 56 ± 1 K. 

Figure [23] presents a plot of versus H/M (an Arrott plot) isotherms near Tc. The 
isotherm that most closely goes though the origin is the one closest to Tc, giving for this 
case a value of 58 K, consistent with the results of the Cp and p(T) measurements. Figure [2^ 
shows magnetization versus external field data along 3 different crystallographic directions: 
[100], [110] and [HI], at 2 K. All of these data sets are consistent with a low temperature FM 
ground state with moderate anisotropy. The spontaneous longitudinal magnetic moment in 
zero applied external field, estimated by the extrapolation of the magnetization curves back 



23 



3 
LL 



1.5 



1.0 



0.5 



55 K ^3 ^ TbFe^Zn^^ 



57 K 

■ 

58 K 

■ 

■ 

. 59 K 



0.0 L.!^ 





. ■ 60 K 



1 

H/M (kOe/|a„ 



FIG. 23: Arrott plot of magnetic isothermals for TbFe2Zn2o- 
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FIG. 24: (Color online) Field dependent magnetization of TbFe2Zn2o along three principle axes at 
2 K. The three lines represent the calculated results based on molecular field approximation are 
all clustered near 9 and appear as a single line. The dashed lines and the values present the 
extrapolate of the magnetization curves and the estimated spontaneous magnetic moments along 
three directions. 



to if = 0, yield M([110]) = S.O/i^, M([lll]) = 6.6/iB, and M([100]) = 5.7^^^. The ratio 
of them is very close to 1 : a/2/3 : a/1/2. Such behavior indicates that the spontaneous 
magnetic moments along [111] and [100] directions can be understood as the projection of 
the one along the easy axis, [110]. At 2 K, the saturated moment at 55 kOe along the easy 
axis, [110], is S.l/i^, O.O/i^ less than the value associated with the Hund's ground state. 
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FIG. 25: (a) Temperature dependent M/H for DyFe2Zn2o {H = 1000 Oe); (b) Cp] (c) p and 
dp/dT. Upper inset: magnetic part of Cp. Lower inset: magnetic entropy. 

2. DyFe2Zn2o 

The low field thermodynamic and transport properties of DyFe2Zn2o are shown in Fig. [251 
The temperature dependent magnetization data (Fig. [25] a) suggest a FM transition below 
50 K. The specific heat data show a kink associated with magnetic ordering (Fig. [2^ b). 
which can be seen more clearly after the subtraction of the non-magnetic background (upper 
inset) and indicates Tc = 45 ± 1 K. This FM transition temperature is further confirmed by 
a weak change in slope in p{T) (associated with the low temperature loss of spin disorder 
scattering), indicating Tq = 45 ± 2 K. Given that the loss of spin disorder scattering in 
intermetallics often scales with de Gennes parameter {32!], the feature we find in p(T) below 
Tc becomes fainter and fainter as R progresses from Gd to Tm. These values of Tq are 
consistent with the result of the Arrott plot analysis, from which a value of Tc = 45 ± 1 K 
can be inferred(Fig. [26]) . 

It is worth noticing that the specific heat data show a faint shoulder near 10 K, which 
appears to be a broad peak after the background subtraction, and is coincident with a slope 
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FIG. 26: Arrott plot of magnetic isothermals for DyFe2Zn2o- 
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FIG. 27: (Color online) Field dependent magnetization of DyFe2Zn2o at 2 K along three princi- 
ple axes. The solid lines represent the calculated result based on molecular field approximation 
(see data analysis part below). The dashed lines and the values present the extrapolate of the 
magnetization curves and the estimated spontaneous magnetic moments along 3 directions. 



change feature in p(T) data. As seen below, such anomaly below Tq in Cp and p(T) data 
also appears for the members of R = Ho, Er and Tm. Those anomalies are likely due to the 
magnetic excitation energy spectrum associated with the Hund's rule multiplet of ions 
in their FM states. (Further discussion will be presented below.) 

The 2 K field dependent, magnetization isotherms for DyFe2Zn2o are shown in Fig. 
[271 Compared to TbFe2Zn2o, the magnetization curves for DyFe2Zn2o reveal a slightly 
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FIG. 28: (a) Temperature dependent M/H for HoFe2Zn2o {H = 1000 Oe); (b) C^; (c) p and 
dp/dT. Upper inset: magnetic part of Cp data. Lower inset: magnetic entropy. 



more complicated, anisotropic behavior. The magnetization along [100] direction manifests 
one metamagnetic phase transition near 12 kOe. Above this transition, the magnetization 
along [100] direction is essentially the same as that for the field along the easy [111] axis. 
The spontaneous longitudinal magnetization along the three directions, M([lll]) = 9.1//^, 
M([110]) = lAfiB, and M([100]) = 5.3/xb, have a ratio very close to 1 : : ^flj?>. These 

results indicate that M([110]) and Af ([100]) can be seen as the projection of M([lll]). The 
metamagnetic phase transition along [100] can be understood as the process of a classical 
spin reorientation in a cubic symmetry coordination. As in the case for GdFe2Zn2o and 
TbFe2Zn2o, the saturated moment of DyFe2Zn2o at 55 kOe, 9.5//s, is slightly less than the 
value of the Hund's ground state value, 10//^- 
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FIG. 29: Arrott plot of magnetic isothermals for HoFe2Zn2o- 



3. HoFe2Zn2o 

Figure [28] presents the low field thermodynamic and transport data from measurements 
on HoFe2Zn2o. The anomalies associated with the FM transition in HoFe2Zn2o in the specific 
heat and resistivity data are relatively weak. The specific heat anomaly can be associated 
with Tc ~ 28 K, and the dp/dT data show faint anomaly at this temperature (Fig. [281) . 
The Tc value is determined as 28 ± 1 K from Cp data, as well as 29 ± 1 K from p(T) data. 
This determinate Tq value is consistent with the result of the Arrott plot analysis (Fig. [2^ . 
which gives Tc = 28 ± 1 K. 

The low temperature magnetic isotherms for HoFe2Zn2o (Fig. [30jl manifest similar, but 
obviously larger, anisotropy to the ones for TbFe2Zn2o- The ratio of the spontaneous mag- 
netization, M([110]) : M([lll]) : M([100]) = 9.1^^ : : G.l/is is close to the ratio of 
1 : a/2/3 : a/1/2. This ratio is consistent with the projection of the local moment from the 
easy [110] axis onto the [111] and [100] axes. In the external field of 55 kOe, the magneti- 
zation along the easy axis, [110], reaches the value of 9.9 fiB, very close to the value of the 
Hund's ground state, lO/i^. 
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FIG. 30: (Color online) Field dependent magnetization of HoFe2Zn2o at 2 K along three princi- 
ple axes. The solid lines represent the calculated result based on molecular field approximation 
(see data analysis part below). The dashed lines and the values present the extrapolate of the 
magnetization curves and the estimated spontaneous magnetic moments along 3 directions. 
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FIG. 31: (a) Temperature dependent M/H for ErFe2Zn2o {H = 1000 Oe); (b) Cp] (c) p and dp/dT. 
Upper inset: magnetic part of Cp data. Lower inset: magnetic entropy. 
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FIG. 32: Arrott plot of magnetic isothermals for ErFe2Zn2o- 
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The low field thermodynamic and transport properties of ErFe2Zn2o are shown in Fig. 
|3T1 The specific heat data show a kink near 18 K (Fig. [31] b), which can be seen more 
clearly after the background subtraction (upper inset) and indicates Tc = 18 ± 1 K. The 
resistivity data show no clear anomaly at this temperature. (Fig. [3T|) The released magnetic 
entropy reaches 21 J/mol K at Tc, 90% of the one associated with the Hund's ground state 
of Er^"*", i?lnl6 (Fig. lower inset). Although p(T) data manifest no anomaly at Tc, we 
will see below that the weak anomaly associated with magnetic ordering can be blow up 
after the background [p(T) for LuFe2Zn2o] subtraction. The Arrott plot for ErFe2Zn2o (Fig. 
I32l) . although showing non-linear, isothermal curves, demonstrates Tc = 17 ± 0.5 K with 
little ambiguity. The non-linear feature is not unexpected for the 4/ local moment systems 



associated with the CEF induced anisotropy. 
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The magnetic anisotropy of ErFe2Zn2o is reminiscent of that of DyFe2Zn2o: both have the 
same easy and hard magnetization orientations, [111] and [110] respectively, as well as the 
metamagnetic transition along the [100] direction (Fig. [331) . The ratio of the spontaneous 
longitudinal magnetic moments, M([lll]) : M([110]) : M([100]) = 7.4//^ : 5.9/xb : 4.2/iB is 
also close to the ratio of 1 : a/2/3 : a/1/3. 
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FIG. 33: (Color online) Field dependent magnetization of ErFe2Zn2o at 2 K along three princi- 
ple axes. The solid lines represent the calculated result based on molecular field approximation 
(see data analysis part below). The dashed lines and the values present the extrapolate of the 
magnetization curves and the estimated spontaneous magnetic moments along 3 directions. 

5. TmFe2 Zn2o 

The low field magnetization, specific heat and resistivity data for TmFe2Zn2o are shown 
in Fig. [311 The temperature dependent magnetization data suggest a FM transition below 
10 K (Fig. [31] a). However, the specific heat data for TmFe2Zn2o only manifest one broad 
peak at 4.5 K (Fig. [3l|b), which is less like the anomalies associated with Tq for R = Gd - 
Er, and more like a Schottky anomaly associated with a CEF splitting. The resistivity data 
also show anomaly below 5 K (Fig. [311 c). However, at this point, it is difficult to determine 
whether this anomaly is associated with the magnetic ordering or the CEF splitting of the 
4/ electrons of Tm^+ ions. As we can see below, after the subtraction of the nonmagnetic 
background, the anomaly associated with the loss of the spin disorder scattering can be seen 
more clearly. 

For TmFe2Zn2o, the Arrott plot analysis provides the reliable criterion for Tq determi- 
nation. Figure [35] shows that Tq can be determined as 5.5 ± 0.5 K without any ambiguity. 
At this temperature, the magnetic entropy is 15 J/mol K, only 70% of the value of fully 
released entropy of Hund's ground state of Tm^+, _Rlnl3 (Fig. [311 upper inset). 

The low temperature magnetic isotherms for TmFe2Zn2o manifest the same easy and 
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FIG. 34: (a) Temperature dependent M/H for TmFe2Zn2o {H = 1000 Oe); (b) Cp- (c) p and 
dp/dT. Upper inset: magnetic part of Cp data. Lower inset: magnetic entropy. 



hard axis as Tb and Ho members, [110] and [111], respectively (Fig. 1361) . The spontaneous 
longitudinal magnetic moments along the three principle axes are all close to 4/^^- Such a 
result may be due to the relatively low value of Tc, which makes the spontaneous magnetic 
moment less anisotropic at 2 K. The saturated moment along the easy axis reaches 6.2/is 
at 55 kOe, O.S/i^ less than the value of the Hund's ground state, T/i^. 



D. YbFe2Zn2o and YbCo2Zn2o 



The physical properties of YbFe2Zn2o and YbCo2Zn2o have been discussed in ref. jll j. 
They manifest typical heavy Fermion behavior with large electronic specific heat, 520 
mJ/mol and 7900 mJ/mol K^, respectively. Figure EZI shows temperature dependent 
susceptibility and resistivity data for these two Yb compounds. The susceptibility data for 
YbFe2Zn2o manifest a broad, Kondo-type peak about 20 K, indicating a clear loss of local 
moment behavior, whereas the susceptibility for YbCo2Zn2o shows CW behavior down to 
1.8 K (Fig. [12]), associated with the effective moment value /Xe// = 4.5/xb- Above ~ 50 K, 
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FIG. 35: Arrott plot of magnetic isothermals for TmFe2Zn2o- 
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FIG. 36: (Color online) Field dependent magnetization of TmFe2Zn2o at 2 K along three princi- 
ple axes. The solid lines represent the calculated result based on molecular field approximation 
(see data analysis part below). The dashed lines and the values present the extrapolate of the 
magnetization curves and the estimated spontaneous magnetic moments along 3 directions. 



x{T) for YbFe2Zn2o manifests a CW behavior with an effective moment of 4.7fiB, close to 
the value of the Hund's ground state of Yb^+, 4.5/2^ (see Fig. [2T1) . The resistivity data 
for YbFe2Zn2o show a broad shoulder about 30 K, whereas for YbCo2Zn2o, the resistivity 
data shows a Kondo resistance minimum about 50 K and a clear coherent peak about 2 K. 
These apparently different behaviors for these two Yb-based heavy fermion compounds with 
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FIG. 37: Temperature dependent M/H (a) and resistivity (b) for YbFe2Zn2o and YbCo2Zn2o 
{H = 1000 Oe). 

virtually identical structures has been explained as the result of significantly different Kondo 
temperatures: Tk = 33 K and 1.5 K for Fe and Co compounds, respectively. 

IV. DATA ANALYSIS AND DISCUSSION 

As shown in Fig. [3H] a, the Tq values of RFe2Zn2o compounds (R = Gd - Tm) scale 
fairly well with the de Gennes factor, dG = {gj — 1)^J(J + 1), consistent with a RKKY 
interaction. All of the dc values for each compounds are smaller than their respective Tc 
values (for R = Er and Tm, the values of 9c are even negative). These small 9c values are 
consistent with the deviation of x{T) from the CW law (Fig. [2T1) . As observed in the case 
of pseudo-ternary compounds Gd2:Yi_3;Fe2Zn2o, such deviation can be explained as a result 
of increasing coupling between the local moments embedded in the strongly temperature 
dependent, polarizable matrix, YFe2Zn2o or LuFe2Zn2o- [idI 




Previous studies show that the magnetization of GdFe2Zn2o at base temperature is nearly 
isotropic with a slightly reduced saturated moment (~ O.bfis less than the value of the 
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FIG. 38: (Color online) Tq and 9c (a), the maximum value on Ap (b) with respect to the de 
Gennes factor for RFe2Zn2o (R = Gd - Tm). 



Hund's rule ground state of Gd^+). For R = Tb - Tm, the magnetization anisotropy at 
base temperature is significant, and correlates with the easy and hard axes of the respective 
RCo2Zn2o analogues. Such behavior indicates the anisotropy of the RFe2Zn2o (R = Tb - Tm) 
compounds may mainly be due to the CEF effect on the R^"*" ions. The M{H) curves at 2 K 
manifest divided behavior with R = Tb, Ho and Tm on one hand and R = Dy and Er on the 
other hand: for R = Tb, Ho and Tm, the magnetization processes are gradual along all 3 
principal axes; for R = Dy and Er, the magnetization data along [100] direction shows clear 
and sharp metamagnetic transition. Both types of magnetization processes (gradual increase 
and metamagnetic transition) are common for the FM ordered 4/ local moments with CEF 
anisotropy associated with the R in a cubic point symmetry, and can be understood in terms 
of the purification of the CEF split 4/ electronic wave function due to the Zeeman effect of 
the external field, and the rotation of the local moment, Given that Tb^"*" and Tm^+, as 
well as Dy^"*" and Er^"*" ions have same total 4/ electronic Hund's rule ground state quantum 
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number (J = 6 and 15/2 respectively), the similar magnetic anisotropy indicates similar 
CEF effect for the two sets of rare earth ions, respectively. 

In order to better understand the magnetic anisotropy of RFe2Zn2o compounds (R = 
Tb - Tm), the CEF effect acting on the R ions must to be considered. However, multiple 
difficulties associated with the strongly polarizable back ground [Y(Lu)Fe2Zn2o] as well as 
the strong magnetic interaction, make the determination of the CEF parameters hard. For 
example, in order to reduce the magnetic interaction, the magnetic R^+ ions were placed 
into a dilute coordination, Ra;Yi_^Fe2Zn2o or R^Lui_a;Fe2Zn2o- A FM ground state has been 
found even for very dilute magnetic R concentration: it was found that Tbo.o5Yo.95Fe2Zn2o, 
Dyo.o5Yo.95Fe2Zn2o as well as Hoo.iYo.9Fe2Zn2o manifest FM ordering above 2 K. For such 
small X, the background subtraction (magnetization and/or specific heat of YFe2Zn2o or 
LuFe2Zn2o), as well as the uncertainty of x, make the fitting process unreliable. 

On the other hand, due to the very similar R coordination and the lattice parameters for 
Fe and Co series, the CEF parameters determined from RCo2Zn2o compounds should be close 
to those for the RFe2Zn2o compounds, for respective R members. Figure [39] shows that the 
anisotropic magnetization isotherms for the pseudo-ternary Fe compounds, Ero.iYo.9Fe2Zn2o 
and Tmo.iYo.9Fe2Zn2o, which still manifest a paramagnetic state at the base temperature, 
are close to the calculated results from the determined CEF parameters of the respective 
Co compounds. The calculated results also fairly well mimic the crossing behavior of the 
magnetization along [110] and [100] directions for R = Er, as well as along the [111] and 
[100] directions for R = Tm. For all three directions, the calculated results are slightly 
larger than the experimental ones, which is most likely due to the ±0.02 uncertainty of the 
nominal x value. The larger magnetization for Ero.iYo.9Fe2Zn2o than the calculated results 
below 10 kOe is consistent with residual FM interactions between the Er^"*" local moments. 

The magnetization along the three axes for the all Fe compounds were calculated based on 
the molecular field approximation in a self-consistent manner. In the single-ion Hamiltonian 
for the R^+ ions (Eqn. [3]), with the molecular field approximation, the magnetic interaction 
term is written as: 



exc 



(6) 



where Hm is the molecular field. It obeys the self-consistent condition: 




(7) 
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FIG. 39: (Color online) Field dependent magnetization for Ero.iYo.iFe2Zn2o (a) and 
Tmo.iYo,iFe2Zn2o (b) along three principle axes at 1.85 K. The solid lines present the calculated 
magnetization by using the single-ion Hamiltonian and the CEF coefficients determined from the 
respective Co members. 



/ f\ ^ 'ZoJnexp{-E„/kBT) 

where J„ and En are the eigenvalues and eigenenergies of the nth eigenfunction; A is the 
molecular field constant which can be obtained from the ordering temperature: A = 

Me// 

The calculated magnetizations were compared with the experimental results in Figs. EH 
[271 [30} [331 and [361 All these calculated magnetization values are obviously larger than the 
experimental results. This difference is probably due to (i) the molecular field approximation 
over-estimating the molecular field constant as well as the internal field, and (ii) the induced 
moments from the Fe site aligning in an antiparallel manner with respect to the R^"*" local 
moments (as in the case of GdFe2Zn2o pil]). 

Figure HOl shows the magnetic part of specific heat as a function of T /Tq for RFe2Zn2o (R 
= Gd - Tm). The magnetic ordering temperature {Tq) of R = Gd - Er members manifests 
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FIG. 40: (Color online) Magnetic part of specific heat versus T/Tq for RFe2Zn2o (R = Gd - Tm). 

itself as the position of maximum slope, with a decreasing sharpness as R varies from Gd 
to Er. TmFe2Zn2o does not appear to have any anomaly in the Cm data at Tq. Below 
Tc, the data sets for R = Dy - Tm show a broad peak, which shifts closer to its Tq as R 
varies from Dy to Tm, whereas the data for GdFe2Zn2o show no explicit peak. If the broad 
peaks are corresponding to the magnetic excitation energy spectrum associated with CEF 
effect, then the relative positions of these peaks to Tc, to some extend, indicate the ratio of 
the energy scales of the CEF splitting (for a single ion) to the magnetic interaction. The 
shift of the peak position as R varies from Dy to Tm indicates that the energy scale of the 
magnetic order decreases relative to the CEF splitting. Such phenomena is consistent with 
the analysis on the magnetic part of entropy: as shown in the insets of Figs. [221 ESI EHl [31] 
and [311 Tb, Dy and Ho compounds manifest fully released Sm at their Tc; whereas Er and 
Tm compounds still release part of Sm above their Tc, which indicates that, unlike R = Gd 
- Ho members, the CEF splitting for the 4/ electronic configuration of the Tm^+ and Er^+ 
may extend above magnetic ordering temperature. 

Based on the assumption that the Fe and Co series have similar CEF splitting (for a 
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FIG. 41: (Color online) Single-ion CEF splitting energy levels for RCo2Zn2o (R = Tb - Tm). The 
arrows present the Tq values for RFe2Zn2o with respective R. 



similar R ions), the comparison between the magnetic ordering temperature and the CEF 
splitting for different R ions is qualitatively diagrammatized in Fig. HU The levels represent 
the single ion, CEF splitting of the Hund's ground state of 4/ electronic configuration of 
R^"*", determined from RCo2Zn2o and the arrows represent the Tc values of RFe2Zn2o. The 
Tq value is comparable with the highest energy level of CEF splitting for R = Ho. For R 
= Er and Tm, the Tq values is about | and | of the highest CEF levels, respectively. This 
diagram, though it cannot be used to determine the precise energy splitting of the RFe2Zn2o 
compounds (the CEF levels have been strongly modulated and mixed by the interaction 
energy), is qualitatively consistent with the specific heat measurements, and indicates that, 
at least for TmFe2Zn2o, the CEF energy splitting already happens well above its Tq- In 
summary, is appears plausible that, due to extremely similar liganal environments, equivalent 
members of the RFe2Zn2o and RCo2Zn2o series have similar CEF splitting schemes. 

Further insight can also be gained from a careful revisiting of the transport data. The 
total resistivity of RFe2Zn2o (R = Gd - Yb) can be written as: 

p{T) = Po + Pph (T) + p.mag (T) , (9) 

where Pmag is associated with scattering from the 4/ moments and the spin fluctuation of 
itinerant electrons. As seen in Fig. |12] a, for the whole series above 250 K, the resistivity 
data sets show essentially linear behavior with slopes differing by less than 12%, within 
the estimated dimension error (±10%) of these bar-like-shape samples. These similar, high 
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FIG. 42: (Color online) (a): p versus T, (b): Ap versus T, (c): Ap versus T/Tq for RFe2Zn2o (R 
= Gd - Tm). Inset: the blow up Ap data for TmFe2Zn2o- The arrow presents the FM ordering 
temperature. 



temperature behaviors indicate that, in the high temperature hmit, the magnetic scattering 
is saturated, and the phonon scattering is essentially similar for the whole series (due to the 
very dilute nature of the R ions). Therefore, the magnetic contribution to the resistivity can 
be estimated by (1) subtracting residue resistivity, po (2) normalizing the high temperature 
slope of all p{T) to that of LuFe2Zn2o and then (3) subtracting the piu{T) — Pluo data from 
the normalized data. The result is written as: 

dp| 

mathrmR I 

Ap(r) = (PR - PK0)^^P^ - (PLu - PLuo). (10) 

~dT^ I275K 

As shown before, the subtraction background pi„(T) already includes the scattering 
associated with the spin fluctuation of itinerant electrons. Thus, Ap will not only include 
the scattering from the 4/ moments, but will also include scattering associated with the 
interaction between the 4/ moments and itinerant electrons. Figure |12] b and c show Ap 
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versus temperature, as well as normalized temperature (T/Tc) for R = Gd - Tm. For R = 
Gd - Er, a pronounced upward cusp, whose height decreases from Gd to Er, is centered about 
Tq, whereas TmFe2Zn2o manifests a broad feature and only very weak anomaly around its 
Tc (see the blow-up inset of Fig. W2\i . As shown in Fig. [3H] b, the maximum values of 
the cusps for different R scale with the de Gennes factor, which indicates that the decrease 
of Ap with T below Tc is the result of a loss of spin disorder scattering of conduction 
electrons, associated with the 4/ local moment. However, as found in the pseudo-ternary 
compounds Gda;Yi_j.Fe2Zn2o, 12] the decreasing behavior of Ap with increasing T above Tc 
is more conspicuous and must come from a different conduction electron scatterin g p rocess 



(simple models of p(T) due to magnetic scattering cannot explain this anomaly |35l . |36| ) . 
Giving that RT2Zn2o compounds only manifest this behavior when the local moments are 
embedded in the highly polarizable background (GdCo2Zn2o does not show this behavior 
[lo|). this anomaly is thought to be associated with the spin fluctuation of the 3d electrons. 
Also appearing in the resistivity of RC02 (R = Gd - Tm)[37t], this decreasing behavior of 
Ap with increasing T above Tc has been explained as the result of the increase of the 
spin fluctuation of 3d electrons, which is provided by the increasing, nonuniform fluctuating 
4:f-d electron exchange interaction, as the temperature approaches Tc in the paramagnetic 
state. Since both Y(Lu)Co2 and Y(Lu)Fe2Zn2o are classical examples of NFFLs, such an 
anomaly could be associated with these strongly correlated electron systems. On the other 
hand, considering that the Hund's ground state of Tm'^^ has been significantly split above 
FM ordering for TmFe2Zn2o, it is not unexpected that the conduction electron scattering 
process manifests a different behavior associated with the CEF effect. 

Finally, the nearly FM compounds: YFe2Zn2o and LuFe2Zn2o also merit some further 
discussion. Shown in Fig. [3l the low field susceptibility (H = 10 kOe) manifests a maximum 
near 6 K and 8 K for YFe2Zn2o and LuFe2Zn2o respectively. Such a maximum in the tem- 
perature dependent susceptibility also appears for other examples of nearly FM compounds. 
For example, Pd manifests Tmax ~ 70 K 20 1; YC02 and LuCo2 manifests Tmax ~ 100 Kj2ll|; 



and TiBe2 manifests Tmax ~ 10 K 221]. Another interesting phenomena in nearly FM mate- 



rials is the so-called itinerant electron metamagnetism (lEM), which is an applied magnetic 
field induced, first order, phase transition between a paramagnetic state and spin polar- 



ized state 



381 ]. Experimentally, lEM has been observed for YC02 and LuCo2 around 70 T. 



39l . 1401 ] Within the framework of Landau theory, the maximum in temperature dependent 
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susceptibility is thought to be related to lEM. 41j The magnetic part of the free energy AF 
can be writen as the function of the magnetic moment M: 

AF = -aM^ + -bM^ + -cM^ , (11) 
2 4 6 ^ ^ 

where a, b and c are the Landau expansion coefficients. 

As shown by Shimizu 411], the condition for the existence of JEM is: a > 0, 6 < 0, 
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c > and ^ ^ 20' Within the framework of the spin fluctuation theory, Yamada 
generalized this work by introducing a temperature dependent function of the mean square 
amplitude of spin fluctuations. These theoretical works demonstrated that the existence of 
lEM is associated with the maximum in x{T) by means of the factor of |f , which can be 
estimated as: 

Furthermore, the lEM can only happen below Tmax- These results seem to be consistent 
with the experimental results in various itinerant electronic systems. 43] 

According to the Eqn. [121 the values of ^ can be estimated as 310 and 72 for YFe2Zn2o 
and LuFe2Zn2o respectively {M/H ~ x{T) at 10 kOe), which are much larger than the region 
of the existence of lEM, indicating that lEM may not exist. Indeed, recent measurements 
on a part of the LuFe2Zn2o sample used for the magnetization data in Fig. [3] in a pulse 
magnetic field up to 55 T at 0.3 K, show no evidence of metamagnetic transition. In nearly 
FM materials, no evidence of lEM appears for TiBe2, 4J], which also manifests a relative 
low value of T^ax- From these points of view, Y(Lu)Fe2Zn2o and TiBe2 may represent the 
examples of NFFLs different from YC02 and LuCo2. 

This lack of an lEM sheds further light on the magnetic properties of the local moment 
bearing, RFe2Zn2o (R = Gd - Tm) compounds. As shown before, all the members manifest 
2nd order paramagnetic to ferromagnetic phase transitions. This behavior is different from 
that seen in the RC02 (R = Gd - Tm) system: the magnetic phase transitions of R = Dy 
- Tm members for RC02 are 1st order whereas R = Gd and Tb members have 2nd order 
transitions sl. This difference is not difficult to explain in Landau theory: unlike Y(Lu)Co2, 
the host of Y(Lu)Fe2Zn2o lack of ability to show lEM and therefore can not be induced to 
show metamagnetic transition by any molecular field associated with the 4/ local moments. 
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V. SUMMERY 



RFe2Zn2o and RCo2Zn2o (R = Gd - Lu, Y) demonstrate diverse magnetic properties. 
The conspicuous differences between these two related series are mainly associated with the 
conduction electron polarizability of the host (non-magnetic) compounds. YFe2Zn2o and 
LuFe2Zn2o manifest similar, nearly ferromagnetic properties. When the 4/ local moments 
are embedded in this highly polarizable medium, the RFe2Zn2o (R = Gd - Tm) series shows 
highly enhanced FM ordering temperatures. In contrast, YCo2Zn2o and LuCo2Zn2o manifest 
normal, Pauli paramagnetic behaviors. In a related manner, GdCo2Zn2o and TbCo2Zn2o 
show low temperature AFM ordering, and the magnetic properties for RCo2Zn2o (R = Dy - 
Tm) are more strongly influenced by the CEF splitting of the R ions (the dominant energy 
scale). The CEF coefficients determined for the Co series are consistent with the observed 
anisotropics of the Fe series, indicating that the CEF splitting of the R-ions is similar for 
the all Zn Frank-Kasper polyhedra that fully encompass the R-site. On the other hand, 
YbFe2Zn2o and YbCo2Zn2o manifest different heavy Fermion behaviors. 
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APPENDIX A 

Figure US] shows the magnetization (at H = 1000 Oe) and zero applied field resistivity 
for three batches of TbFe2Zn2o, which were synthesized from different initial molar ratio of 
starting elements, Tb:Fe:Zn = 2:3:95, 2:4:94 and 2:5:93. The ferromagnetic ordering tem- 
peratures, determined as 52 ± 2 K, 56 ± 1 K and 67 ± 2 K by resistivity measurements 
for the three samples, increase as the initial concentration of Fe in the solution increases. 
We also found similar features for R = Gd and Er, but the variation of their Tq values 
are less than in the Tb case (Fig. H4l) . Comparative, single crystal x-ray diffraction mea- 
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FIG. 43: (Color online) (a): Temperature dependent M/H for TbFe2Zn2o {H = 1000 Oe) from 
different initial growth molar ratio of starting elements; (b) temperature dependent p in zero 
applied field. 
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FIG. 44: Tq values for different growth of RFe2Zn2o (R = Gd - Tm) samples with respect to the 
de Gennes factor. 
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surements performed on the samples, albeit inconclusive, indicated that the crystallographic 
difference were perhaps due to subtle variations of occupancy of the Fe site.[l8| The main 
difficulty associated with the x-ray diffraction measurement was completely resolving the 
mixed site occupancies for Zn and Fe which have similar atomic number values. Recently, 
carefully prepared, pieces of TbFe2Zn2o samples with same geometric dimension, from the 
starting elements, Tb:Fe:Zn = 2:3:95 and 2:5:93, were measured using single-crystal neu- 



tron diffraction. 19j This measurement showed that the Fe site has ~ 1% deficiency for the 
2:3:95 sample. All these crystallographic measurements indicate the sensitivity of the mag- 
netic properties to the small disorder for RFe2Zn2o in this family which has local moments 
submerged in a highly polarizable conduction electron matrix. 

APPENDIX B 

As shown in Fig. [1], the distance between the rare earth ion and the Zn NNs, as well as 
NNNs, is close to 3 A; whereas the distance with the next next nearest neighbors (NNNN, 
6 Zn in 48/ site) is larger than 5 A. Due to this isolated, cage-like coordinate of rare earth 
ions, the effect of ions other than this CN-16 Frank-Kasper polyhedron can be neglected in 
the calculation of the CEF coefficients, based on the point charge model. 

Appreciably, the ligands of the rare earth ions in the C-15 Laves compounds (RNi2) form 
a same polyhedron, whose CEF coefficients have been calculated by B. Bleaney j45|, based 
on the point charge model. Therefore, one can directly cite the results: 

3m_7e^V^,^ (Bl) 
5S = 4f-S^-^V'-'>(^IMM>. (B2) 



16 V 363i?J 81i?2 / 
where Zic and ^26 are the charge of the NN and NNN ions {Zi = Z2 = 2 for Zn^^), Ri 
and i?2 is the distance between the R ion and the two sets of ions, (r^) and (r^) are the 
mean fourth and sixth powers of the electronic radius for the 4/-electrons, and (3 and 



are the Steven multiplicative factors. Extracting the values of (r^) and (r^) from ref. 46l |. 



j3 and 7 values from ref. 31|, and Ri and R2 values from the results of single crystal X-ray 
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diffraction 




one can calculate the and Bq values (shown in Table [IT|) . 
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